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I. INTRODUCTION 

The invariance under the Lorentz symmetry group is a 
cornerstone of modern physics, and is strongly supported 
by observations. In fact, all the experiments carried out 
so far are consistent with it PQ, and no evidence to show 
that such a symmetry must be broken at certain energy 
scales, although it is arguable that such constraints in 
the gravitational sector are much weaker than those in 
the matter sector [2|. 

Nevertheless, there are various reasons to construct 
gravitational theories with broken Lorentz invariance 
(LI) [3i. In particular, our understanding of space-times 
at Plank scale is still highly limited, and the renomaliz- 
ability and unitarity of gravity often lead to the violation 
of LI 0j. One concrete example is the Horava theory of 
quantum gravity [5], in which the LI is broken via the 
anisotropic scaling between time and space in the ultra¬ 
violet (UV), 

t^b~ z t, x' L ^b~ 1 x l , (* = 1,2,3), (1.1) 

where z denotes the dynamical critical exponent. This 
is a reminiscent of Lifshitz scalars in condensed matter 
physics [5], hence the theory is often referred to as the 
Horava-Lifshitz (HL) quantum gravity at a Lifshitz fixed 
point. The anisotropic scaling Cl 3 provides a crucial 
mechanism: The gravitational action can be constructed 
in such a way that only higher-dimensional spatial (but 
not time) derivative operators are included, so that the 
UV behavior of the theory is dramatically improved. In 
particular, for z > 3 it becomes power-counting renor- 
malizable Bl% The exclusion of high-dimensional time 
derivative operators, on the other hand, prevents the 
ghost instability, whereby the unitarity of the theory is 
assured [5]. In the infrared (IR) the lower dimensional 


*The corresponding author 

tElectronic address: Anzhong'WangObaylor.edu 


operators take over, and a healthy low-energy limit is pre¬ 
sumably resulted 1 . It is remarkable to note that, despite 
of the stringent observational constraints of the violation 
of the LI [I], the nonrelativistic general covariant HL 
gravity constructed in El is consistent with all the solar 
system tests mm and cosmology [14j Qh]. In addition, 
it has been recently embedded in string theory via the 
nonrelativistic AdS/CFT correspondence [TB]. Another 
version of the HL gravity, the health extension m , is 
also self-consistent and passes all the solar system, astro- 
physical and cosmological tests 2 . 

Another example that violates LI is the Einstein-aether 
theory, in which the breaking is realized by a timelike vec¬ 
tor field, while the gravitational action is still generally 
covariant [T8]. This theory is consistent with all the solar 
system tests ^TS] and binary pulsar observations 1511 . 

However, once the LI is broken, speeds of particles can 
be greater than that of light. In particular, the dispersion 
relation generically becomes nonlinear [14] . 

E 2 = ° 2pp2 ( x+ai iw) + a 2 iw) ) ’ (L2) 

where E and p are the energy and momentum of the 
particle considered, and c p , ct* are coefficients, depend¬ 
ing on the species of the particle, while M* denotes the 
suppression energy scale of the higher-dimensional oper¬ 
ators. Then, one can see that both phase and group ve- 


1 It should be emphasized that, the breaking of LI can have sig¬ 
nificant effects on the low-energy physics through the interac¬ 
tions between gravity and matter, no matter how high the scale 
of symmetry breaking is [9]. Recently, Pospelov and Tamarit 
proposed a mechanism of SUSY breaking by coupling a Lorentz- 
invariant supersymmetric matter sector to non-supersymmetric 
gravitational interactions with Lifshitz scaling, and showed that 
it can lead to a consistent HL gravity m- 

2 In fact, in the IR the theory can be identified with the 
hypersurface-orthogonal Einstein-aether theory m in a particu¬ 
lar gauge [3] 1191120| , whereby the consistence of the theory with 
observations can be deduced. 
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locities of the particles are unbounded with the increase 
of energy. This suggests that black holes may not exist 
at all in theories with broken LI, and makes such theo¬ 
ries questionable, as observations strongly indicate that 
black holes exist in our universe (22 ■ 

Lately, a potential breakthrough was the discovery 
that there still exist absolute causal boundaries, the so- 
called universal horizons , in the theories with broken LI 
[23lj . Particles even with infinitely large velocities would 
just move around on these boundaries and cannot es¬ 
cape to infinity pMj . The universal horizon radiates like 
a blackbody at a fixed temperature, and obeys the first 
law of black hole mechanics [251 ■ 

Recently, we studied the existence of universal hori¬ 
zons in the three well-known black hole solutions, 
the Schwarzschild, Schwarzschild anti-de Sitter, and 
Reissner-Nordstrom, and found that in all of them uni¬ 
versal horizons always exist inside their Killing horizons 
[251J . In particular, the peeling-off behavior of the globally 
timelike khronon field u ^ was found only at the universal 
horizons, whereby the surface gravity K pee iing is calcu¬ 
lated and found equal to E3, 

kuh = - u a D a (u\( x ) , (1.3) 

where and D^ denote, respectively, the time transla¬ 
tion Killing field and covariant derivative with respect to 
the given space-time metric g^ (fi, v = 0,1, 2, 3). For the 
Schwarzschild solution, the universal horizon and surface 
gravity are given, respectively, by (26], 

^Sch. _ sch. _ f 2V /2 }_ (1 4 x 

n UH - 4 > ^peeling 1 3 y r ’ 

where r s denotes the Schwarzschild radius. 

In this paper, we shall study the formation of the uni¬ 
versal horizons from realistic gravitational collapse of a 
spherically symmetric star 3 . To be more concrete, we 
shall consider such a collapsing object in the Einstein- 
aether theory [T5] , To make the problem tractable, we 
further assume that the effects of the aether are negligi¬ 
ble, so the space-time outside of the star is still described 
by the Schwarzschild solution, while inside the star we 
assume that the distribution of the matter is homoge¬ 
neous and isotropic, so the internal space-time is that of 
the Friedmann-Robertson-Walker (FRW). Although the 
model is very ideal, it is sufficient to serve our current 
purposes, that is, to show explicitly that universal hori¬ 
zons can be formed from realistic gravitational collapse. 

Specifically, the paper is organized as follows: In Sec. 
II, we shall present a brief review on the definition of 
the stationary universal horizons, and then generalize it 
to dynamical spacetimes. This is realized by replacing 


3 Here “realistic” means that the collapsing object satisfies at least 
the weak energy condition EH). 


Killing horizons by apparent horizons [251 [50j , and in the 
stationary limit, the latter reduces to the former. In Sec. 

III, we study a collapsing spherically symmetric star with 
a finite thickness in the framework of the Einstein-aether 
theory. When the effects of the aether are negligible, the 
vacuum space-time outside the star is uniquely described 
by the Schwarzschild solution and the junction condition 
across the surface of the star reduces to those of Israel 
m ■ Once this is done, we find that the khronon equa¬ 
tion can be solved analytically when the speed of the 
khronon is infinitely large, for which the sound horizon 
of the khronon coincides with the universal horizon. It is 
remarkable that this is also the case for the Schwarzschild 
solution [26] , for which the universal horizon and surface 
gravity are given by Eq.( |1.4| ). The paper is ended in Sec. 

IV, in which our main conclusions are presented. 


II. DYNAMICAL UNIVERSAL HORIZONS AND 
BLACK HOLES 


A necessary condition for the existence of a universal 
horizon of a given space-time is the existence of a globally 
time-like foliation B H31 EH- This foliation is usually 
characterized by a scalar field cf>, dubbed khronon TS\ . 
and the normal vector u^ of the foliation is always time- 
like, 

u x u x = -1, (2.1) 

where 

^ = X = ~9 a& ^d^. ( 2 . 2 ) 

In this paper, we choose the signature of the metric as 
(—1,1,1,1). It is important to note that such a defined 
khronon field is unique only up to the following gauge 
transformation, 

4> = X W, (2-3) 

where iF((f>) is a monotonically increasing (or decreasing) 
and otherwise arbitrary function of </>. Clearly, under the 
above gauge transformations, we have = u^. 

The khronon field <j> is described by the action [18| . 


V 0 (D^u v f + V 0 ( D ^Z) 2 


+ c 3 (D^u 1 ') (D^u^) - c A a^a h 


(2.4) 


where Cj’ s are arbitrary constants, and a^ = u a D a u 
The operator D M denotes the covariant derivative with re¬ 
spect to the background metric g as mentioned above. 
Note that the above action is the most general one in the 
sense that the resulting differential equations in terms of 
u/j, are second-order m- However, when u^ is written in 
the form of Eq.(2.2), the relation 


n^i/Daiip] 0 , 


(2.5) 
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is identically satisfied. Then, it can be shown that only 
three of the four coupling constants Cj are independent. 
In fact, from Eq.( |2.5[ ) we find [151. 

+ ( D a up ) (. D a u /3 ) - {D a up) (D p u a ) = 0 . 

(2-6) 

Then, one can always add the term, 

A S<f, = c 0 j \J\g\ d 4 xA (2.7) 

into Stj,, where Cq is an arbitrary constant. This is effec¬ 
tively to shift the coupling constants Cj to c', where 

c'l = Cl + Co, C 2 = C2, C3 = C3 — Co, C4 = C4 — Co- (2.8) 

Thus, by properly choosing co, one can always set one of 
Cj (i = 1, 3,4) to zero. However, in the following we shall 
leave this possibility open. 

The variation of S$ with respect to cj) yields the 
khronon equation, 


DpA* = 0, 


(2.9) 


where m, 


A " = 
= 

J a „ = 


(SZ+U^Uv) ^ 

yx * ’ 

J 7 " + c 4 a 7 D 1 ' it 7 , 

(W + Vb5“^ + C3^ 

- c i u a u p g^)Dpu''. ( 2 . 10 ) 


A. Universal Horizons in Stationary and 
Asymptotically flat Spacetimes 

In stationary and asymptotically flat spacetimes, there 
always exists a time translation Killing vector, £**, which 
is timelike asymptotically, 

C A Ca < 0 , (2.11) 

for r —> 00 . A Killing horizon is defined as the existence 
of a hypersurface on which the time translation Killing 
vector becomes null, 


C A Ca = 0 . (2.12) 

On the other hand, a universal horizon is defined as the 
existence of a hypersurface on which becomes orthog¬ 
onal to Ufj,, 


«aC a = 0 . 


(2.13) 


Since is timelike globally, Eq.(2.13) is possible only 


when ^ becomes spacelike. This can happen only inside 
the apparent horizons, because only in that region 
becomes spacelike. 


B. Universal Horizons in Non-Stationary 
Spacetimes 


To study the formation of universal horizons from grav¬ 
itational collapse, we need first to generalize the above 
definition of the universal horizons to non-stationary 
spacetimes. For the sake of simplicity, in the rest of 
this paper we shall restrict ourselves only to spherical 
space-times, and its generalization to other spacetimes is 
straightforward. 

The metric for a specially symmetric space-time can 
be cost in the form, 

ds 2 = gijdx'dx^ + R 2 (x l ) dfl 2 , ( i,j = 0,1), (2.14) 

in the spherical coordinates, = (a; 0 , x 1 ,9, tu) , (fi = 
0,1, 2,3), where dfl 2 = d0 2 + sin 2 ddzu 2 . 

The normal vector n M to the hypersurface R = Cq is 
given by, 

n r = = * 2*0 + ( 2 - 15 ) 

where C 0 is a constant and i?j = dR/dx l . Setting 

= %Ri - S 4 R 0 , (2.16) 

we can see that ^ is always orthogonal to n M , 

C X nx = 0. (2.17) 


For spacetimes that are asymptotically flat there al¬ 
ways exists a region, say, R > R^, in which 
and are, respectively, space- and time-like, that is, 
n ^n^\ R>R > 0 and C m CmI,r>,r < 0- An apparent hori¬ 
zon may form at Rah, at which n M becomes null, 

nxn li?=i? AH = ® 5 (2.18) 


where Rah < Roo- Then, in the internal region R < 
Rah, becomes timelike. Therefore, we have 


( > 0, R> Rah, 
n x n x = < = 0, R = Rah, 
l < 0, R < Rah- 


Since Eq. (2.17) always holds, we must have 


( < 0, R> Rah, 
CaC A = s = 0, R = Rah, 
l > 0, R < Rah, 


(2.19) 


( 2 . 20 ) 


that is, ^ becomes null on the apparent horizon, and 
spacelike inside it. 

We define a dynamical universal horizon as the hyper¬ 
surface at which 


ux< * \r=R uh 


( 2 . 21 ) 


Since u^ is globally timelike, Eq. (2.21) is possible only 


when is spacelike. Clearly, this is possible only inside 
the apparent horizons, that is, Rjjh < Rah- 

In the static case, the apparent horizons defined above 
reduce to the Killing horizons, and the dynamical univer¬ 
sal horizons defined by Eq. (2.21) are identical to those 
given by Eq.(2.13). 
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III. GRAVITATIONAL COLLAPSE AND 
FORMATION OF UNIVERSAL HORIZONS 


Let us consider a collapsing star with a finite radius 
Ry (r), where r denotes the proper time of the surface of 
the star. To our current purpose, we simply assume that 
the space-time inside the star is described by the FRW 
flat metric, 


ds 2 _ = —dt 2 + a 2 (t) (dr 2 + r 2 c?0 2 ) , (R<Ry(t)), (3.1) 

where R = a(t)r is the geometric radius inside the col¬ 
lapsing star. From Eq. (2.15), the normal vector to 
the hypersurface R = Co takes the form, 


n li =d(t)r5°+a(t)6l, (3.2) 

where a = da/dt. Then, the corresponding vector ^ 
reads 


of the observers that are comoving with the collapsing 
surface of the star. In the current case, we have t = t. 
Then, we find 


f?s(r) = a(r )rs, 

(i - v 2 - 2 R-sv -1 = 0, (3.8) 


where a dot denotes the derivative with respect to r. 
However, since in the present case we have r = t, so we 
shall not distinguish it from that with respect to t, used 
above. The extrinsic curvature tensor on the two sides of 
the surface defined by, 


Ka b = 


- d 2 x± 

.d£ a d£ b 


+ r 


0 /> 


dx± dx± - 

~d£?~dt? J’ 


(3.9) 


has the following non-vanishing components [32] 4 


C = a(t)$°u - a{t)r5l. 


(3.3) 


According to Eq.(2.18) the apparent horizon locates at, 

1 


tah = - - 




(3.4) 


Note that in the collapsing case we have a < 0. 

The spacetime outside the collapsing star is vacuum. 
In the Einstein-aether theory pa, if we consider the case 
where the effects of the aether field is negligible, then 
the vacuum space-time will be that of the Schwarzschild, 
and in the ingoing Eddington-Finkelstein coordinates 
(v, R , 9, vj) the metric takes the form, 


_l_ v Mil 
TT = ~b~R^’ 

I<g g = sin~ 2 9K~ lp = a(r)rY, 


(3.10) 


K+ = sin- 2 0K+ v = {Rz-2M)v-RzR- Si , 

where are the normal vectors defined in the two faces 
of the surface 13.61) , 


= 8(r - r s ) 

' vrv r\ U Qi 5 


dx c _ 

d(R-Ry(v)) 

dx°i 


= e R _ dRxjv) , v 
dv a 


(3.11) 


2 M 
~R 


ds 2 + = — ( 1 — ‘CC- ) dv 2 +2dvdR+Fi 2 dfl 2 , (R > R^(t)), 


(3.5) 

where M denotes the total mass of the collapsing system, 
including that of the star surface, and has the dimension 
of length L, that is, [ M ] = L. The surface E of the star 
can be parameterized as, 


r — ty, = 0 (in V ), 

R-Ry(v) = 0 (in V + ), (3.6) 

where is a constant, when we choose the internal co¬ 
ordinates (t, r, 0 , w) are comoving with the fluid of the 
collapsing star. On the surface of the collapsing star, the 
interior and exterior metrics reduce to 


From the Israel junction conditions E2, 


[Aafc] dab [A] — 87TT a fo, 


(3.12) 


we can get the surface energy-momentum tensor r a g, 
where [ K ab ]~ = K+ b - K~ _[K}^ = g ab [K ab ]~ , and g ab 
can be read off from Eq.(3.7), where a,b = r, 6,m. In¬ 
serting Eq.( |3.10 1 into the above equation, we find that 
Tab can be written in the form 


r ab = aw a w b + r](9 a 0 b + zu a ^b), (3.13) 

where w a , 9 a and w a are unit vectors defined on the sur¬ 
face of the star, given, respectively, by w a = 6^, 0 a = 
R^ e a , vo a = Ry sin 95™, and 


ds 2 _ | 


r=rj: 


-dt 2 + a 2 (t)r|dU 2 
ds 2 + \ R=Rs ( v ) 


/ 2 M 2 dRY(v) \ 

\ Rz(v) dv ) 

+Ry(v) 2 d£l 2 
-dr 2 + R s (t)W, 


dv 2 


(3.7) 


where Ry(t) is the geometric radius of the collapsing 
star, v is a function of r, where r denotes the proper time 


a 

V 


1 a 

4irrYd 4ira 
v v 

8irrYd 8 ttv 


Mil 

2ir r^a 2 
1 

87T TyCL 


ii 

47TTYa ’ 

d Mil 
87 ra 87 t r|a 


(3.14) 


2 ’ 


4 Note the sign difference of the first term of K+t between the 
one obtained here and that obtained in m, because of t he s ign 
difference of the cross term dvdR in the external metric (3.5). 


























5 


here a is the surface energy density of the collapsing star, 
and rj its tangential pressure. Physically, they are of¬ 
ten required to satisfy certain energy conditions, such as 
weak, strong and dominant [28], although in cosmology 
none of them seems necessarily to be satisfied [ 88 ] . 

On the other hand, inside the collapsing star, the 
khronon can be parametrized as, 


vA = \Jl + a 2 V 2 Sg + VSi, 

= — \/l + o, 2 V 2 Sq + a 2 V5i, (3.15) 


where V = V(t, r ) is determined by the khronon equation 
(2.9), which now reduces to 


where 


A r r + 


2 A r 


+A + 


3 a{€)A t 

a(t) 


= 0 , 


(3.16) 


A‘ = Ci23-4i + CuA* 2 + c 13-83 

A r = ' /l \f r V2 A\ A e = A* = 0 , (3.17) 

arv 


with Cab = C a + ct, c a bc = c a + c b + c c , and 

A{ = 9 , ^ ( 2 a 4 P 5 (-l + 2 r 2 ad) 

r-{ 1 + a z V z ) 

+r{rV" + V'{2 + ra 2 \J 1 + a 2 V 2 V)) 

+V (—2 — 3r 2 d 2 + 3 r 2 aa + 5 r 2 a\J 1 + a 2 V 2 aV' 
+r 2 a 4 V 2 + 2ra 2 Vl + a 2 V 2 {V + rV')) 
+a 2 P 3 (—4 + r{—2ra 2 + 4 ra\J 1 + a 2 V 2 aV' 
+ci(7ra + 2a\/l + a 2 V 2 {V + rV')))) 
+ra 4 V 4 (2V' + r(V" + a{ 5dV + aV))) 
+ra 2 V 2 (4V' + r(2V" + a(7dl> + aV)))), (3.18) 

Al = - —^——(4rVa 2 +ra 5 V 4 (2Va + 5aV) 

r(l + a z V z ) 

+2a 3 V 2 {2V 2 Vl + a 2 V 2 a + 2rV(a 
+ V 7 ! +a 2 V 2 aV') + 5raV) + a(4V 2 y /1 + a 2 V 2 cl 
+ rV (2d + 5\/1 + a 2 V 2 aV') + 5rdV) 

+ra 6 V 4 V + a 4 V 2 ( 4rV 3 d 2 + V 2 {2V' + rV") 
+2V\J\ + a 2 V 2 {V + rV') + 2 rV) 

+a 2 (8rV 3 d 2 + V 2 (2V’ + rV") + V{rV 2 
+2^1 + a 2 V 2 (V + rV')) 

+r(\/l + a 2 V 2 V'V + V))), (3.19) 

Al = 2P 2 (1 + a 2 V 2 )(a 2 — ad). (3.20) 


coupling constants Cj. However, when C 14 = 0, we obtain 
a particular solution, 

V(t,r) = -770 a(t) = a 0 e~ Ht , (3.21) 

a(t) 

where Vq, H and ao are integration constants with 
H > 0, ao > 0. It is remarkable to note that C 14 = 0 cor¬ 
responds to the case in which the speed of the khronon 
becomes infinitely large = C 123 /C 14 —> 00 , a case that 
was also studied in [8 [23] [26]. 

From the definition of the dynamical universal horizon 
Eq.(2.21) and considering Eqs.(3.3), (3.15) and (3.21), we 


find that the collapse always forms a universal horizons 
inside the collapsing star, and its location is given by, 




^V 4 + 4V 2 a 2 H 2 e~ 2Ht - V 2 ' 


(3.22) 


From Eq.(|3.8|), on the other hand, we find that, 

(3.23) 
back into 


v = 


—r|aa + 

2 M — arY. 


Subs tituti ng it together with a{r) = aoe H 
Eqs. (3.14), we obtain 




V = 


87ri?s 


M-Ry, - 2H 2 R 3 : 

RyQ 


- 1 


where 


Ry = ao^se 


-Hr 


Q = 




(3.24) 


(3.25) 


Obviously, to have both a and 77 real, we must assume 
that Ry > R^ m , where R^f ln is a root of the equation, 


1 + H 2 Rl - 


2 M 

Ry: 


= 0 . 


(3.26) 


When the star collapses to the point Ry ( TMin ) = -Rvf ln , 
the tangential pressure diverges, whereby a space-time 
singularity (with a finite radius) is developed. This rep¬ 
resents the end of the collapse, as the space-time beyond 
this moment is not extendable. 

It can be shown that the weak energy condition a > 
0 , a + r] > 0 can always be satisfied by properly choosing 
the free parameters involved in the solution, before the 
formation of the universal horizon. In particular, from 


Eq.(3.24) we evan see that er is always non-negative, and 

1 


cr + r] = 


87 tRyQ 


g-(^~ 1 


(3.27) 


Thus, for Ry > 3 M, the weak energy condition is always 
satisfied. When Ry < 3 M, it is also satisfied, provided 
that Q > 3 M/Ry — 1, or equivalently 


Here a prime denotes the derivative with respect to r. 
It is found very difficult to solve Eq.(3.4) for any given 


h 2 rI 


4 M 
Ry: 


> 


9 M 2 


(3.28) 
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Clearly, by properly choosing the free parameters in¬ 
volved in the solution, this condition can hold until the 
moment where the whole collapsing star is inside the uni¬ 
versal horizon. However, at the end r = TMin of the 
collapse this c onditi on is neces sarily violated, as can be 
seen from Eqs.(3.26) and (3.281. Using the geometric ra¬ 
dius R = ra(t) = ra 0 e~ Hi inside the collap sing star, we 
find that the apparent horizon given by Eq.(3.4) and the 


universal horizon given by Eq.(3.221 can be expressed as, 


Ruh = ruria{t) 


— a^e 


— Ht 


Rah = rAHa(t) = 


y/VJ 1 + 4 V 0 2 a 2 0 H 2 e~ 2Ht - U 0 2 
1 


H 


(3.29) 


In Fig. |T] we show one of the cases, in which the free 
parameters are chosen as = 1 ,M = l,ao = 3, Vo = 
0.6, H = 1.5. Then, we find that the weak energy condi¬ 
tion holds until the moment t = 0.624554, at which we 
have Ruh = 1.24387 > R% = 1.1756. That is, the weak 
energy condition holds all the way down to the moment 
when the whole star collapses inside the universal hori¬ 
zon, as it is illustrated clearly in Figj2| In this figure, 
three horizontal lines, R = 2M, 3M/2, R^ m are also 
plotted, where Ruh = 3M/2 is the universal horizon in 
the Schwarzschild space-time [cf. Eq. (1.4)]. For the cur¬ 
rent choice of the free parameters, the universal horizon 
is not continuous. In fact, when the star collapses to the 
moment t = t a , where t a is given by i? s (f 0 ) = 3M/2, 
the universal horizon jumps from Ruh (to) to 3M/2, as 
shown more clearly in Fig. [3] Physically, this is be¬ 
cause that the surface shell of the collapsing star has 
non-zero mass. The collapse ends at t = t^ ln , where 
Rs(t^ Im ) = at which the surface pressure of the 

collapsing star becomes infinitely large. It is remark¬ 
able that R^ m generically is different from zero, that 
is, the collapse generically forms a space-time singular¬ 
ity that has finite radius. The corresponding Penrose 
diagram is given in Fig. |4j In this figure the location 
of the event horizon denoted by the straight line EH is 
also marked, although it can be penetrated by particles 
with sufficiently large velocities, and propagate to infin¬ 
ity, even they are initially trapped inside it. However, 
this is no longer the case when across the universal hori¬ 
zon. As explained above, once they are trapped inside 
the universal horizon, they cannot penetrate it and prop¬ 
agate to infinities, even they are moving with infinitely 
large velocities. As a result, an absolutely black region 
is (classically) formed from the gravitational collapse of 
a massive star, and this region is black, even in theories 
that allow instantaneous propagations! 


IV. CONCLUSIONS 

In this paper, we have first generalized the definition 
of a stationary universal horizon to a dynamical one, 



FIG. 1: The solid (black) line represents the surface energy 
density a, while the dashed (red) line denotes the quantity 
o + r). When plotting these curves, we set rs = 1, M = 
l,ao = 3, Vo = 0.6,7/ = 1.5. At the moment t = 0.624554, 
we have o + r) = 0, R s = 1.1756, and Ruh = 1.24387. 



FIG. 2: Universal horizon Ruh, the surface radius Rs of the 
collapsing star vs t for rn = 1, M — 1, oo — 3, Vo = 0.6, H = 
1.5. The horizontal lines R = 2 M and R = 3JU/2 are, respec¬ 
tively, the Killing and universal horizons of the Schwarzschild 
vacuum solution. The collapse ends at the moment when 
Rs = A(f m , at which point the pressure of the surface of the 
collapsing star becomes infinitely large, whereby a space-time 
singularity with finite thickness is developed. 


by simply replacing Killing horizons by apparent ones. 
Then, we have constructed an analytical model that rep¬ 
resents the gravitational collapse of a spherical symmet¬ 
ric star with finite thickness, and shown explicitly that 
dynamical universal horizons can be formed from such a 
“realistic” gravitational collapse. Here “realistic” is re¬ 
ferred to as a gravitational collapse of a star with a finite 
thickness that satisfies at least the weak energy condition 

m- 

To have the problem tractable, we have assumed that 
the star consists of an anisotropic and homogeneous per¬ 
fect fluid and that outside the star the space-time is 
vacuum in the framework of the Einstein-aether the¬ 
ory [IS]. When the effects of the aether field is negli¬ 
gible, the vacuum space-time is uniquely described by 
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FIG. 3: Gravitational collapse of a spherically symmetric star 
with its radius R^{t), which divides the space-time into two 
regions, V ± . The curved line Ruh denotes the universal hori¬ 
zon formed inside the star, while the vertical straight line R = 
3M/2 is the universal horizon of the Schwarzschild vacuum 
solution. When the star collapses inside the Schwarzschild 
universal horizon R = 3M/2, the universal horizon suddenly 
jumps from RuH(t 0 ) to 3M/2, because of the non-zero mass 
of the collapsing surface of the star. 


the Schwarzschild solution [18]. Even in this case, solv¬ 


ing the khronon equation (2.9) inside the star is still very 
complicated. Instead, we have further assumed that the 
velocity of the the khronon is infinitely large, so the sound 
horizon of the khronon coincides with the universal hori¬ 
zon. Then, we have found that an analytical solution 
exists for the star made of the de Sitter universe, and 
shown explicitly how a dynamical universal horizon can 
be formed, as can be seen from Figs. i-S 

Although such a constructed model serves our current 
purpose very well, that is, to show that universal hori¬ 
zons can be indeed formed from realistic gravitational 
collapse, it would be very interesting to consider cases 
without (some of) the above assumptions, specially the 
case in which the space-time outside of the star is not 
vacuum, so that the star may radiate, when it is collaps¬ 
ing. 

In addition, although in this paper we have consid¬ 
ered the formation of the universal horizons only in the 
framework of the Einstein-aether theory, it is expected 
that our main conclusions should be true in other theo¬ 
ries of gravity with broken LI, including the HL gravity 

DU 113- 
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